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Abstract 

We present the first experimental limits on high-g^ contributions to charmless semileptonic B 
decays of the form expected from the Weak Annihilation (WA) decay mechanism. Such contribu- 
tions could bias determinations of \Vub\ from inclusive measurements oi B ^ X^lv. Using a wide 
range of models based on available theoretical input we set a limit of rwA/r6-»w < 7.4% (90% 
confidence level) on the WA fraction, and assess the impact on previous inclusive determinations 
of \Vub\. 
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A precise determination of the magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element Vub with well- understood uncertainties is one of the highest priorities in 
heavy- flavor physics. Recently, significant experimental progress in the determination of \ Vub\ 
through inclusive measurements of semileptonic B decays has been achieved with pioneering 
work such as the incorporation of 6 — 57 spectral information for improved modelling [l| 
and the use of large -B-tagged samples 0, 1^ • Theoretical advances include the quantitative 
evaluation of the leading and subleadingcontributions to the partial B — >■ X^ii^ width in 
restricted regions of phase space 0, IE HiS 1^ • Several concerns remain, one of which is the 
the "Weak Annihilation" (WA) contribution [1, HO, III to the total b uiu rate. 



The WA contribution arises from a four-quark operator at order {A/rrib)^ in the heavy 
quark Operator Product Expansion (OPE) in which, crudely, the b quark annihilates within 
the B meson. Due to an enhancement of IGvr^ relative to other (9(m^^) operators, WA rate 
estimates, l^, 11| 

rwr.-.«o.o3(^)^(^). (1) 

are a sizeable few percent. A non-zero WA contribution requires violation of QCD factoriza- 
tion, parameterized above by B2 — Bi. Little is known about the scale of this violation since 
it is fundamentally non-perturbative. Because WA is expected to be concentrated in phase 
space near k, ml, its relative importance is magnified by kinematic requirements that 
accentuate the high region to isolate b ^ u from the large b ^ c background. Limiting 
WA is thus important for understandin g th e precision of inclusive \Vub\ determinations. 
In this Letter, we describe a search 2| for a b ^ uiu rate at large g^, such as would 



be expected from WA. We use the 15.5 fb ^ of data collected at the T(45') resonance with 
the CLEO II |i3], CLEO II. V Q and CLEO III detectors at the Cornell Electron 



Storage Ring (CESR). The analysis employs the missing momentum and missing energy 
techniques used in analyses of semileptonic moments [la. [l7l | and originally developed in 
studies of exclusive charmless semileptonic decay H, 1^ 20|- The procedure takes advantage 



of the near hermeticity of the CLEO detectors to estimate the four momentum of a single 
missing neutrino from the missing four momentum Pmiss = (-Emiss,Pmiss) in an event, where 

-Emiss = 2i?beam — Ei and ^miss = PCM — Y. Pi- 
All detector configurations provide acceptance over more than 90% of the full An solid 
angle for both charged particles (momentum resolution of 0.6% at 2 GeV/c) and photons 
(typical 7r° mass resolution of 6 MeV). Charged particles are assigned the most probable mass 
based on a combination of specific ionization measurements, either time-of-fiight or Cerenkov 
radiation angle measurements, and the relative spectra of vr^, and p from B decay. To 
optimize the missing-momentum resolution, the selection is optimized to identify a single 
track for each charged particle originating from the B and B decays with high efficiency, 
and to measure energy deposited in the Csl calorimeter that is not associated with charged 
tracks or their interaction products. The track subset choice is based on event topology 
rather than individual track quality since high efficiency with minimal double counting is 
more important than use of only well-measured tracks. Energy clusters located within 8 cm 
of tracks determined to project into the Csl calorimeter, or consistent with being "split off" 
from a matched shower, are excluded from the energy and momentum sums. 

Electrons satisfying p > 400 MeV/c are identifed over 90% of Att using the ratio of cluster 
energy to track momentum in conjunction with specific ionization (dE/dx) measurements 
in the main drift chamber. Time-of-flight or Cerenkov measurements provide additional 
e^/K^ separation in the momentum range with ambiguous dE/dx information. Particles 
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in the polar angle range |cos^| < 0.85 that register hits in counters beyond 5 interaction 
lengths are accepted as signal muons. Those with | cos^| < 0.71 and hits between 3 and 
5 interaction lengths are used in the multiple-lepton veto discussed below. We restrict 
signal electrons and muons to the range 1.5 < p < 3.0 GeV/c. Within these fiducial and 
momentum regions the selection efficiency exceeds 90%. The probability of misidentifying a 
charged hadron as an electron (muon) is about 0.1% (1%). 

In events with multiple undetected particles, Pmiss represents Pi, poorly and causes re- 
constructed variables in i? — ^ Xciiy decays to smear beyond their kinematic limits into the 
regions of sensitivity for the much rarer B Xud-v process (including WA). This mechanism 
provides the dominant background contribution in this analysis. Therefore, we reject events 
with multiple identified leptons, which are usually accompanied by multiple neutrinos. We 
also reject events with a net charge not equal to zero, or where Pmiss is consistent with 
particles lost down the beam pipe (| cos^missl > 0.9). 

The missing mass squared M^j^^ = i^^j^g— |pmissp, which should be zero (within resolution) 
if only a sole neutrino is missing, provides further background suppression. Because the 
resolution on -Emiss is about 60% larger than that on IpmissL the M^- resolution aj.j2 

miss 

2EraissO'E,T,iss- Requiring M^^^j2E^iss < 0.2 GeV provides a zero mass requirement at roughly 
constant -Emiss resolution that enhances the signal relative to the background by over a factor 
of two. We then take pi, = {Ey,py) = (|pmiss|, Pmiss) in other kinematic calculations. In 
particular, we can calculate g^, the square of the hadronic momentum transfer in semileptonic 
decay, via = {pe+p^Y- The core resolution on is about 0.6 GeV^, with a broad high-side 
tail from events with more than one undetected particle. 

Finally, we ensure that the event is consistent with e~^e^ — > BB decay. We begin with 
the standard CLEO hadronic sample, defined by events with at least six primary tracks 
and a visible energy of at least 20% of the center of mass energy. We suppress continuum 
e~^e~ — >• qq backgrounds and r^r~ backgrounds using the ratio (i?2 = H2/H0) of the second 



to zeroth Fox- Wolfram moments [21| and a sphericity-like variable [2^ that is sensitive to 
momentum fiow perpendicular to the lepton, which is small for continuum events. Of the 
events that satisfy all other criteria and > 2 GeV^, our continuum suppression rejects 
72% of continuum while retaining over 90% of semileptonic B decays. 

We search for evidence of WA or other sources of high-g^ B —>■ Xuiu decays by fitting our 
measured g^ spectra in three lepton-momentum bins. The lowest bin (1.5 < pe < 2.0 GeV/c) 
is dominated by i? ^ X^fz/ and serves to normalize this background. The total B X^ii^ 
yield is determined by the middle (2.0 < pi < 2.2 GeV/c) and highest {p£ > 2.2 GeV/c) bins, 
while the greatest sensitivity for WA-like processes is provided by the highest momentum 
bin. Our g^ spectra for the full sample and for the highest momentum bin are shown in Fig.^ 
Corrections have been applied for continuum background (using data collected just below 
BB threshold) and for events in which hadrons were misidentified as signal leptons (using 
data obtained with a lepton veto folded with measured misidentification probabilities). 

The corrected g^ spectra are fitted with B Xciu and B Xuii^ components obtained 
from simulation. The B Xciv simulation incorporates all available 5-decay data, includ- 
ing measured semileptonic decay form factors. The B X^^iu simulation is based on (i) a 
hybrid model that combines the HQET-based approach of DeFazio and Neubert [24] with 
known exclusive resonances, and (ii) a simple model for WA that refiects both the kinematics 
implied by c/FwA/c^g^ ~ ^(g^ — ^1) ^"^^ the intuitive picture where the "valence" quarks in 
the B meson annihilate and a soft non-perturbative hadronic system Xu materializes. In 
this formulation, the lepton-neutrino pair carries most of the energy (g2 ~ M|), while the 
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FIG. 1: The continuum- and fake-lepton-subtracted spectra (points) for pi > 1.5 GeV/c and 
Pl > 2.2 GeV/c (inset) with components B Xciiy (hght grey), B Xuli' (hatch) and WA 
(with {Mx) = 0.293 GeV, dark grey). The dashed envelope results from systematic variation of 
B XJu. 



TABLE I: Summary of "impact ratios" for some WA models considered. 



xo (GeV) A (GeV) 


{Mx) (GeV) 


-Rtotal (%) 


-Rcndpoint (%) 


Rq'^M^ (%) 


Rm^ (%) 


0.30 


0.01 


0.293 


1.73 ±0.68 ±0.60 


8.24 ±3.04 ±2.34 


3.89 ±1.50 ±1.18 


2.24 ±0.88 ±0.67 


0.30 


0.05 


0.328 


1.40 ±0.69 ±0.58 


6.73 ± 3.11 ± 2.36 


3.17±1.52±1.18 


1.82 ±0.89 ±0.66 


0.30 


0.20 


0.476 


1.36 ±0.89 ±0.82 


5.61 ± 3.52 ± 3.13 


3.06 ±1.97 ±1.77 


1.76 ±1.15 ±1.00 


0.30 


0.30 


0.574 


1.49 ± 1.25 ± 1.16 


5.19 ±4.16 ±3.82 


3.33 ±2.72 ±2.50 


1.91 ± 1.58 ± 1.43 


0.30 


0.50 


0.773 


2.95 ± 1.89 ± 1.89 


7.23 ±4.43 ±4.27 


6.02 ±3.74 ±3.64 


3.51 ±2.24 ±2.13 


0.40 


0.01 


0.342 


1.17 ±0.65 ±0.55 


5.62 ± 2.98 ± 2.28 


2.64 ±1.45 ±1.13 


1.52 ±0.84 ±0.63 


0.40 


0.05 


0.369 


1.07 ±0.67 ±0.58 


5.08 ± 3.06 ± 2.51 


2.42 ±1.50 ±1.24 


1.38 ±0.87 ±0.70 


0.40 


0.20 


0.498 


1.22 ±1.00 ±0.96 


4.85 ± 3.82 ± 3.64 


2.76 ±2.22 ±2.11 


1.58 ± 1.29 ± 1.21 


0.40 


0.30 


0.593 


1.88 ±1.27 ±1.30 


6.18 ±3.99 ±4.01 


4.16 ±2.74 ±2.76 


2.39 ± 1.61 ± 1.58 


0.50 


0.01 


0.392 


0.93 ±0.71 ±0.61 


4.35 ±3.20 ±2.70 


2.10 ±1.58 ±1.34 


1.20 ±0.91 ±0.76 


0.50 


0.05 


0.416 


0.95 ±0.76 ±0.75 


4.36 ±3.34 ±3.44 


2.17 ±1.69 ±1.72 


1.24 ±0.98 ±0.99 


0.50 


0.10 


0.452 


1.03 ±0.84 ±0.94 


4.44 ±3.48 ±4.11 


2.34 ± 1.87 ±2.16 


1.34 ± 1.08 ± 1.26 


0.50 


0.20 


0.534 


1.44 ± 1.07 ± 1.07 


5.27 ±3.76 ±3.73 


3.25 ±2.37 ±2.34 


1.86 ±1.38 ±1.34 


0.50 


0.30 


0.621 


1.58 ±1.36 ±1.32 


4.90 ±4.07 ±3.88 


3.51 ±2.96 ±2.83 


2.01 ± 1.72 ± 1.62 


0.50 


0.50 


0.806 


3.42 ± 2.32 ± 2.35 


7.56 ±4.91 ±4.79 


6.85 ±4.49 ±4.42 


4.04 ±2.72 ±2.64 


0.60 


0.01 


0.442 


0.82 ±0.80 ±0.81 


3.64 ±3.45 ±3.55 


1.87 ± 1.81 ± 1.84 


1.07 ± 1.04 ± 1.06 


0.60 


0.05 


0.465 


1.01 ±0.87 ±0.91 


4.28 ± 3.58 ± 3.84 


2.29 ± 1.96 ±2.07 


1.31 ± 1.13 ± 1.19 


0.60 


0.30 


0.660 


2.67 ± 1.64 ± 1.62 


7.41 ±4.34 ±4.10 


5.82 ±3.47 ±3.32 


3.38 ±2.06 ±1.92 


0.60 


0.50 


0.836 


4.27 ±2.61 ±2.75 


8.71 ±5.07 ±5.10 


8.44 ±4.93 ±5.02 


5.00 ±3.03 ±3.03 


0.75 


0.01 


0.518 


1.32 ±1.04 ±1.10 


5.01 ±3.79 ±4.01 


2.99 ±2.31 ±2.43 


1.72 ±1.34 ±1.40 


0.75 


0.20 


0.641 


3.33 ± 1.79 ± 1.89 


9.20 ±4.65 ±4.70 


7.31 ±3.77 ±3.85 


4.27 ±2.27 ±2.24 


0.75 


0.30 


0.719 


3.43 ±2.09 ±2.09 


8.11 ±4.71 ±4.45 


7.38 ±4.32 ±4.13 


4.31 ±2.61 ±2.43 


0.75 


0.50 


0.886 


4.63 ±2.86 ±3.80 


8.34 ±4.94 ±6.37 


8.97 ±5.28 ±6.85 


5.34 ±3.27 ±4.20 



hadronic system has kinematics at the non-perturbative scale Aqcd- To describe the spectra 
of that soft hadronic system, our implementation introduces a probability density function 
(pdf) that is flat out to a cutoff Xq, where an exponential roll-off with slope A begins. The 
mass Mx and momentum of the hadronic system for a WA decay are drawn independently 
from this pdf, uniquely determining the kinematics. The system is then hadronized into at 
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least two particles or resonances. The kinematics of the ^v pair are calculated assuming the 
V — A structure of the weak current and spin s = for the hadronic system. We examine 
combinations of five xq and six A values, for a total of 30 different WA cases that span a 
wide range of kinematics. 

We perform a separate fit for each WA case, with the b ^ c, b ^ u, and WA rates 
floating independently. The WA rate is not constrained to be positive. Acceptable fits 
to the distributions are obtained for all cases. The B — >■ X^^z/ and non-WA B — 
Xuiiy components dominate, and in no case is the WA component more than two standard 
deviations above zero (combined statistical and systematic). The most significant WA yield 
is obtained for the case with {M^) = 0.293 GeV (shown in Fig. Q), and appears to result 
from an overlap with a downward fluctuation in the distribution of the sample used for 
continuum subtraction. 

From each fit's results for the WA and non-WA b ^ u rates, the ratio R = FwA/rb^u 
is computed for the full phase space (total) and for restricted phase space regions that 
have been used in previous inclusive \Vub\ measurements: > 2.2 GeV/c (endpoint); pi > 
1.0 GeV/c, g2 > 8.0 GeV^ and Mx < 1.7 GeV (g^M^); and > 1.0 GeV/c, Mx < 1.55 GeV 
(Mj). These ratios constrain the extent to which a measured rate can be biased away from 
current theoretical estimates because of a localized WA contribution. The ratios for the full- 
phase-space and endpoint cases are shown in Fig. |21 and the results for a subset of the cases 
considered are summarized in Table HI In each region our results, which are statistics-limited, 
set non-trivial constraints on a localized WA enhancement. 

The primary systematic uncertainties arise from experimental effects related to recon- 
struction of the neutrino, such as the absolute and b ^ c ^ siu rates and spectra, the 
efficiency and resolution for charged particle and photon detection, modeling and rejection of 
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FIG. 2: Fractional size of the WA component for the full phase space (bottom) and restricted to 
Pi > 2.2GeV/c endpoint region (top). The statistical (total) uncertainties are represented by the 
inner (full) error bar. 
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TABLE II: Systematic uncertainties (WA model of Fig. 



Source 


Ai?tot 


Aiitot/i?tot(%) 


7 efficiency 


0.00177 


10.2 


tracking efficiency 


0.00247 


14.3 


resolution 


0.00095 


5.5 


Ptrk resolution 


0.00134 


7.7 


Kl multiplicity 


0.00013 


0.8 


hadronic shower modeling 


0.00118 


6.8 


hadronic shower veto 


0.00065 


3.8 


particle identification 


0.00078 


4.5 


b ^ c ^ s£u 


0.00020 


1.1 


b ^ civ modeling 


0.00349 


20.1 


b ^ uiv modeling 


0.00309 


17.9 


Total 


0.00601 


34.7 



charged hadronic showers, and charged hadron identification 20] . The B Xciu modehng 
systematic estimate includes variations of the branching fractions at levels commensurate 
with recent measurements, and variations of form factors at levels several standard deviations 
from recent average results 23j . The B Xj,£i> modeling systematic includes a variation of 
the inclusive shape function similar to Ref. [l| and variation of the X„ hadronization model. 
In the ratios of the WA component relative to the B XuC-v component, many common 
systematics related to luminosity, fake rates, etc. largely cancel. Table |n] summarizes the 
systematic contributions for the WA model shown in Fig.^ Shifts observed with systematic 
cross checks such as fioating individual components of the b ^ ciu background, fioating in- 
dividual classes of mistakes (e.g., extra Ki or extra z/) in the b ^ ciu background sample, 
and even more extreme variations such as eliminating the Diu and nonresonant b ^ ciu 
components are commensurate with the quoted systematics. 

To limit the bias in rate measurements quantitatively, we parameterize the variation in 
the central value and total uncertainty of each WA fraction (R) as a function of (Mx)- 
90% confidence limits (CL) are then calculated assuming a fiat probability distribution in 
{Mx) over the range we have investigated, resulting in -Rtotai < 7.4%, -Rendpoint < 15.5%, 
Rq2M^ < 14.5% and Rm^ < 8.6%. Limits on a bias of \Vuh\ are half these values. These 
results provide the first concrete constraint on one of the three important uncertainties in 
extraction of \Vub\ for which only more speculative results i^] have existed to date. They 
also place \Vub\ from endpoint analyses on a much stronger footing, where the 8% bound 
(90% CL) we find for a bias in an endpoint-based \Vub\ is much more restrictive than the 



estimate ctwa ~ 10 — 20% j2J|. Studies like these are crucial for inclusive determinations of 
\Vub\ to achieve a 5% precision goal robustly (already achieved statistically). 

In summary, we have obtained the first experimental limits on the potential bias in 
inclusive determinations of \Vub\ from a localized contribution to the distribution, as 
could arise from weak annihilation. The method presented here is one of several, including 
study of semileptonic rates in the charm sector and of the moments of the B semileptonic 
distribution 2^|, that will be needed to understand weak annihilation and its impact upon 
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ub 
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